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Abstract[] Eulophids provide a valuable experimental model system to investigate a wide variety of questions in the 
evolution of parasitic wasps. Notauli[] an important taxonomic character in the family Eulophidae[] is used to 
differentiate subfamilies[] tribes[] genera and species. However] abundant morphological homoplasy in this family has 
made it difficult to confidently identify many eulophid taxa at specific level. Both Cirrospilus ambiguus and Diglyphus 
bimaculatus are considered to have notauli patterns that are intermediate between their respective genera. Determining 
the phylogenetic position of both species may facilitate studying the evolution of notauli pattems. We analyzed CO 1 [] 
1181 and 28S sequences using Bayesian method. C. ambiguus was transferred into the monophyletic Diglyphus[] and 
the monophyly of Cirrospilus[] however[] was not supported. The types of notauli were examined in the context of the 
phylogeny based on 285 gene sequences. The complete notauli in the Eulophidae extending to the hind margin of the 
mesoscutum was inferred to represent the ancestral form. The complete notauli that curves to meet the apex of the 
axillae occurs independently in five clades. Each occurrence represents a unique independent evolution. The 
incomplete form of notauli occurred in four clades[] indicating that the shape independently evolved four times. 

Key words|] Eulophidae[] molecular systematics[] Cirrospilus ambiguus[] Diglyphus bimaculatus|] notauli[] evolution 


behavior]  physiology[] ecological preference 


1 INTRODUCTION 


The Eulophidae] Hymenoptera[] Chalcidoidea[] is 
one of the most speciose families of parasitoid wasps. It 
currently contains 294 genera and 4 288 species 
П Noyes[] 2002[]. Besides attacking eggs and larvae[] 
eulophids can also attack pupae and adults of their 
hosts. These species occur in almost all known 
parasitoid forms[] endoparasitoids and ectoparasitoids|[] 
idiobionts and koinobionts[] solitary and gregarious[] and 
specialists and generalists. With such a diverse 
biology[] eulophids form excellent model systems to 
study the evolution of diverse life styles and numerous 
other questiond] Godfray[] 19941. 

Leafminers[] the hosts of eulophid parasites[] cause 
great damage to crops and flower$] Minkenberg and van 
Lenteren[] 1986[]. One of the important[] effective and 
environmentally friendly methods to control leafminers is 
However[] the incorrect 
identifications of parasitoids and their hosts will hinder 


to use natural enemies. 


their use in biological control. No matter how similar 
morphologically[] the parasitoid species can differ in 





reproductive modes [] Rossler and DeBach[] 19730 
Gordh[] 1982[]. Many species of both genera Cirrospilus 
and Diglyphus of Eulophidae are important parasitoids 
of leafminers[] Zhu et al.[] 2000[] Noyes[] 2002[T] 
including Cirrospilus ambiguus|] Hansson and LaSalle[] 
19961 which attacks Liriomyza trifolii throughout the 
Tropicd] eg[] China[] Indial] Vietnam and Tanzania[] 
[] Noyes[] 2002[] and in South Africa. 

Three notauli shapes occur in the Eulophidae[] 
incompletd] Fig. 11 complete and curving to meet the 
apex of the axilla[] Fig. 2[[] and complete and 
extending to the hind margin of the mesoscutunl] Fig. 
3[]. These patterns are used to differentiate eulophid 
subfamilies[] tribes and generd] Boucek[] 1988[] Schauff 
et аі. 1997[]. Using this character] Peck et al. 

O 1964[] divided the Eulophinae into two subfamilies[] 
Eulophinae and Elachertinae. The tribe Anselmellini[] 
established by Bouéek[] 1988[[] is easily recognized 
with the combination of three characters[] one of which 
is that the mesoscutum has clear-cut[] deep notauli. 
Tribe Ophelimini[] defined by Ashmead] 190410 has 


complete notauli and two short hind tibial spurs. 
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Bouček [| 1988[] proposed the Ophelimini should 
accommodate most genera with deep[] straight notauli. 

Both C. ambiguus and D. bimaculatus have 
similar notaul{] Fig. 20400 which are considered to be 
examples of intermediate morphological forms between 
Cirrospilus and Diglyphus. The complete and curved 
shapes of notauli are one of the key features 
distinguishing both species from all members of the tribe 
Cirrospilinf] Hansson and LaSalle[] 1996[] Zhu et al .[] 
1999[] Zhu and Huang[] 2002[]. Most members of 
Diglyphus can be distinguished from other genera by 
having incomplete notaul{] Zhu et al .[] 20000. 

Notauli patterns are also used to distinguish 
Cirrospilus from Zagrammosoma [] Gordh[] 19780 
Boucek[] 1988[] LaSalle[] 1989[]. In Cirrospilus[] except 
for C. ambiguus|] notauli extend to hind margin of 
mesoscutum[] Fig. 3[]. In Zagrammosomal] once 
treated as a subgenus of Cirrospilus|] Peck et al.[] 
1963[[] notauli are curved and extend to the anterior 
margin of the axillae. Recently[] many authors recognize 
its generic status[] Gordh 1978[] Bouéek[] 1988[] 
LaSalld] 1989[]. Although Cirrospilus ambiguus is 
placed in its genus because of its metallic body color 
and non-vaulted vertex[] in making this decision 
Hansson and LaSalle[] 1996[] acknowledged that the 
inclusion of C. ambiguus might leave Cirrospilus a 
paraphyletic genus. 

It is difficult to identify eulophids based only on 

morphology because of the relatively small number of 
informative characters and high levels of homoplasy 
0 Ubaidillah et а. 2003[]. Fortunately[] modern 
molecular methods offer important alternatives to species 
identification and they can yield insights into the 
evolution of these parasitic wasps. Relatively 
conservative nuclear gene sequences can distinguish 
species and biotypes of parasitoids that are difficult to 
separate morphologically] Hoy et al .[] 20000. Nuclear 
ribosomal 115] can resolve the relationships among 
populations and species[] Avise[] 2000[]. The second 
expansion segment[] D2[] of 28S rDNA is useful for 
resolving relationships among families іп the 
Chalcidoided] Campbell et a/.[] 2000[]. Gauthier et 
al.[] 2000[] found that this region was particularly 
informative at the levels of subfamily and tribe. The 
rapid evolution and selective constraints make 
mitochondrial protein-encoding genes useful molecular 
markers for separating species[] Hebert et al .[] 20030 
Hebert et al .[] 20041. 

The main purposes of this study аге) 10 to 
investigate the generic placement of the two species[] 
C. ambiguus and D. bimaculatus based on а 
phylogenetic evaluation of sequences from three genes[] 
and] 20 to investigate the evolution of notauli patterns 
in the Eulophinae based on 28S sequences. 





Fig. 1 Dorsal view of the thorax 
of Sympiesid] incomplete[] 





Fig. 2 Dorsal view of the thorax of Cirrospilus ambiguus 
П complete and curving to meet the apexof the axilla[] 





Fig. 3 Dorsal view of the thorax of Cirrospilus brevis 
П complete and extending to the hind margin of the mesoscutum[] 


2 MATERIALS AND METHODS 


2.1 Taxa sampling 
All sequenced specimens[] Tables 1[] 2[] were 
identified by the second author based on morphological 
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features. All voucher specimens were deposited in the 
Institute of Zoology[] Chinese Academy of Sciences[] 
Beijind] IZCAS[]. Mymar sp.[] Mymaridae[] was used 
as the initial outgroup to re-root 5096 majority rule 
consensus tree constructed from 285 D2. 
Morphologically] members of Mymaridae have closer 
relationships with other families than Eulophidae in 
Chaleidoidea[] Schauff[] 19840 Gibson[] 1986[]. The 
genus Euderus was used to root the network constructed 
бот CO | and 1151 sequence data. Euderinae are 
considered as the most primitive members of the family 
because they possess eight gastral segments and always 
four or even five funicular segments[] Boucek[] 1988[] 
Gauthier et al .[] 20000. 





Fig. 4 Dorsal view of the thorax of Diglyphus bimaculatus 
П complete and curving to meet the apex of the axilla[] 


2.2 DNA extraction|] PCR[] and sequencing 

Because most specimens were critical point dried 
before DNA extraction[] they were submerged in tri- 
distilled water for 24 h. Genomic DNA was extracted 
using standard phenol-chloroform extraction from a 
single individual[] Sambrook et а!.[] 1989[]. After 
extraction[] the DNA was dissolved in 450 pL TE buffer 
and stored at — 20%. 

A fragment of CO | was amplified with standard 
polymerase chain reaction] PCR[] protocol and directly 
sequenced after purifications. Primers used for CO [ 
PCR and sequencing were as the following] CO | SF 5'- 
TAAGATTTTGATTATT[] АСОСО TA []CC-3'[] COI2191[] 
5'-СССССТААААТТААААТАТАААСТТС-3' CO I 
2613[] 5'-ATTGCAAATACTGCACCTAT-3' and СО | 
2193[] 5'-TTTTTTGGTCATCCAGAAGT-3'. The first 
primer was designed by ourselves[] while others were 
taken from Simon et al.[] 1994[] and Chen et al. 

O 2004 []. Primers used for amplifications and 
sequencing of 1151 were taken from Ji et al.[] 20030 
П 18sF1 and 5p8sBld[]. For amplifying the D2 region[] 
primers were taken from Campbell et al [] 19930 28S 
D2F and D2R[]. All PCRs were carried out in standard 
25 uL reactions. For mitochondrial CO | [] 35 cycles 
were performed. In each cycle[] melting temperature 


was 94^C. for 30 s[] annealing at 48°C for 45 s[] and 
extension at 72°C for 1 min. For nuclear ribosomal 28S 
and ITSI[] PCR used 35 cycles of 94°C. denaturation for 
30 s[] 50% for annealing 30 s and 72% extension for 1 
min[] with an initial denaturation of 95% for 5 min and 
a final extension of 72°C for 7 min. PCR products were 
purified with Minipore protocols [] Qiagen[] and 
sequenced with ABI Big Dye protocols . 

CO [ and 1151 were not always sequenced from 
the same samples[] Table 1[]. The DNA of some 
individuals was exhausted after multiple attempts to 
amplify one of these gene fragments. Thus[] additional 
samples were used to sequence other fragments. 

2.3 Sequence alignment and phylogenetic analyses 

Sequences were initially aligned by Clustal W 

П Thompson et al .[] 1994[] with default parameter$] gap 
opening penalty = 100 gap extension penalty = 5[] delay 
divergent sequences = 40%[] and then manually 
checked in BioEdii[] Hall[] 1999[]. Missing data were 
coded as“ D". 

To select an evolutionary model that best fits the 
data[] hierarchical likelihood ratio tests [] hLRTs[] 

[] Goldman[] 1993[] were implemented by comparing log- 
likelihood scores of 56 models of base substitution using 
Modeltest[] Posada and Crandall[] 1998[]. Bayesian 
inference was used to estimate genealogical relationships 

П Huelsenbeck and Ronquist[] 20010 Buckley et al .[] 
2002[] Nylander et al.[ 2004[] Ronquist[] 2004[]. 
Analyses were conducted using  MrBayes 2. 0 

[] Huelsenbeck and Ronquist[] 2001[]. АП Bayesian 
analyses were initiated with random starting trees. Four 
Markov chains were used and the dataset was run for 2 
x 10Ó generations to allow for adequate time of 

Trees were sampled every 100 

generations. To avoid being trapped on local optimal] 
we ran two independent analyses with different starting 
trees and monitored the fluctuating value of the 
likelihood graphically 0 Huelsenbeck and Bollback[] 
2001[]. We plotted the log-likelihood scores of sample 
points against generation time and determined that 
stationarity was achieved when the log-likelihood values 
of the sample points reached stationarit[] Huelsenbeck 
and Ronquist[] 2001[]. After discarding burn-in 
samples[] the remaining samples were retained for 
generating a 5096 majority rule consensus tree. These 
frequencies of resolving nodes were assumed to reflect 


convergence. 


the true probabilities of the clades under the assumed 
models. Consequently[] Bayesian posterior probabilities 
[] BPP[] of 9596 or greater were considered to be 
significantly supported . 


3 RESULTS 


Nineteen CO | sequences from 19 species and 18 
1181 sequences from 14 species were obtained] Table 
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1. Thirteen 28S D2 sequences of 12 species were also 
obtained [] Table 2[]. АП 285 D2 sequences of 
Eulophidae[] except for those species transferred out of 
the family by Gauthier et а!.[] 200000 were directly 
downloaded from GenBank for analyses. One sequence 
of Cirrospilus]] AF345598[] was discarded for its little 
homology with other eulophids. The aligned sequences 
are available upon request . 

3.1 Sequence composition 


The totals of 745 bp of СО | were easily aligned 
because no indels were detected. Base frequencies 
averaged A = 0.297[] C = 0.120[] G = 0.152[] T = 
0.432. There was a strong adenine and thymind] A-T[] 
bias. Compared with the first and second positions[] the 
average A-T frequency in the third position was 0.9220 
showing very strong A/T bias[] as was typical of insect 
mitochondrial genomes] Simon et al .[] 1994[] Arias and 
Sheppard[] 1996[]. 


Table 1 Samples and localities for parasitoids for which sequences of CO [ and ITS1 were obtained 








s Voucher . Elevation GenBank accession no. 
Species name Locality Sex 
number Оа СО | ITS1 

Aulogmnus 
A. spl TBbsbd 101 f Aulol01 Bangda[] Basu[] Xizang 4 390 Female = DQ390436 
A. sp GXjxsts 2 f Aulo2 Shengtangshan[] Jinxiu[] Guangxi 1 095 Female DQ149154 = 
Cirrospilus 
C. ambiguus VIETNAMhn 105 f Cirr105 Hanoi[] VIETNAM Female DQ390432 DQ390437 
C. diallus TBcd 1 f Cirrl Changdul] Xizang 3 400 Female DQ390434 - 
C. diallus TBcd 102 f Сит102 Changdu[] Xizang 3 400 Female DQ390433 - 
C. diallus QHgem 103 f Сит103 Geermu[] Qinghai 2 880 Female = DQ390438 
C. pictus BJmtg 101 1 f Cirrl01 1 MentougouL] Beijing 1 300 Female = DQ390439 
C. pictus BJmtg 101 2 f Сит101 2 Mentougou[] Beijing 1 300 Female = DQ390440 
C. variagatus SXhyhs 4 m Cirr4 Hengshan[] Hengyuanl] Shanxi 1 475 Male DQ149157 = 
C. vittatus QHgem 104 f Сит104 Geermu[] Qinghai 2 880 Female DQ149158 DQ390441 
Diglyphus 
D. begini TBrkz 147 m Digl47 Rekaze[] Xizang 3 890 Male Е АҮ948107 
D. begini TBls 187 т Dig187 Lasa[] Xizang 3 650 Male DQ149159 - 
D. bimaculatus 'TBjz 128 f 018128 Relong{] Jiangzi[] Xizang 4 700 Female DQ149160 AY948108 
D. bimaculatus TBls 130 f Dig130 Lasa[] Xizang 3 650 Female DQ149161 AY948109 
D. crassinervis QHdl 119 f Digl19 Deling[] Qinghai 2 900 Female DQ149162 AY948110 
D. isaea ОНа! 108 m Digl08 Qilian[] Qinghai 2 790 Male DQ149174 AY948092 
D. ава ОНа! 126 f Digl26 Qilian[] Qinghai 2 790 Female DQ149176 AY948094 
D. minoeus ТВј 101 1 m "16101 Relong{] Jiangzi[] Xizang 4 700 Male DQ149192 AY948111 
D. minoeus TBjzrl 102 f Digl02 Relong{] Jiangzi[] Xizang 4 700 Female = АҮ948112 
D. расһупеигиѕ QHdl 121 т 016121 Delinghal] Qinghai 2 900 Male DQ149193 AY948113 
D. pulchripes QHgem 146 f Digl46 Geermu[] Qinghai 2 880 Female DQ390435 DQ390444 
Eulophus 
E. sp2 TBed 102 f Eulo102 Changdul] Xizang 3 400 Female DQ149197 DQ390443 
Entedon 
E. sp2 BJmtg 1280 Ent2 Mengtougou[] Beijing 1 280 Female DQ149194 - 
E. sp3 SCbx 2200 Ent3 Baoxing ] Sichuan 2 200 Female DQ149195 = 
Euderus 
E. spl BJzhwy 101 f Eud101 Beijing Botanical Garden[] Beijing 100 Female DQ149196 DQ390442 


After manual alignment[] the ITS1 sequences 
consisted of 1 005 bp positions[] which included partial 
18S sequences in 3' and 5.8S sequences in 5' end. 
Base frequencies averaged A = 0.242[] C = 0.25[] G = 
0.276[] T2 0.231. The average frequency of A/T was 
47.3%[ being roughly equal to that of С/С. Aligned 
285 D2 sequence yielded 486 bp. Base frequencies 
averaged A = 0.221[] C = 0.305[] G = 0.139[] T = 
0.336. 

3.2 Phylogenetic analyses 
The TVM + I + G model was chosen as best-fit 
model for CO | dataset by hierarchical likelihood ratio 


test] – InL = 4 644.5327[]. The model allowed for 
equal transition rates and nucleotide substitution with 
invariable sites plus a G distribution of variable-rate 
site. 

After approximately 300 000 generations[] the log- 
likelihood values of each sampled tree stabilized. The 
last 17 001 sample trees were used to estimate the 5096 
majority rule consensus tree and the BPP4] Fig. 50. 
The tree was moderately well resolved. Unresolved 
relationships occurred within the genus Diglyphus only 

O Clade DU. Both Cirrospilus ambiguus and D. 
bimaculatus were resolved in the clade containing other 
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species of Diglyphus. Howeverl] the support for this 
group] Clade D[] was low[] BPP = 54[]. The other 


species of Cirrospilus fall into three clades[] AL] B and C 
[] Fig. 50. 


Table 2 Samples and localities for parasitic wasps for which sequences of 28S were obtained 








Specimen Elevation GenBank accession 
Species name Locality Sex 
no. 0 то по. 
Aulogmnus 
A. spl TBbsbd 101 f Aulol01 Bangdal] Basu[] Xizang 4 390 Female DQ390414 
Cirrospilus 
C. ambiguus VIETNAMhn 105 f Сит105 Напо | VIETNAM Female DQ390415 
C. pictus BJmtg 101 f Сит101 Mentougou[] Beijing 1 300 Female DQ390416 
C. variagatus SXhyhs 4 m Cirr4 Hengshan[] Hengyuan[] Shanxi 1 475 Male DQ390417 
Diglyphus 
D. albinervis QHdlhb 132 f 016132 Delingha[] Qinghai 2 900 Female DQ390418 
D. bimaculatus TBjz 128 f 4700 015128 Relong{] Лапе | Xizang 4 700 Female DQ390419 
D. bimaculatus TBls 130 f 3650 016130 Lasa[] Xizang 3 650 Female DQ390420 
D. inflatus GSwx 131 m 010131 Wenxian[] Gansu 960 Male DQ390421 
D. isaea SCgzkd 174 f Digl74 Kangding{] Сапа] Sichuan 2 400 Female DQ390422 
D. minoeus h6 DigH6 London[] U. K. Female DQ390423 
D. pachyneurus ZJtms 116 f Digl16 Tianmushan[] Zhejiang 350 Female DQ390424 
Eulophus 
E. sp2 TBed 102 f Eulo102 Changdu[] Xizang 3 400 Female DQ390428 
E. sp3 GSzq 103 f Eulo103 Zhouqiu[] Gansu 3 400 Female DQ390429 
Elachertus 
E. fenestratus SC 1 m Elachl Manigange[ ] Dege[ ] Sichuan 3 965 Male DQ390425 
Elasmus 
E. sp ВЈму 5 f Elasl Beijing Botanical Garden[] Beijing 100 Female DQ390427 
E. spl LNdl 102 m Elas102 Dalian[] Liaoning 300 Male DQ390426 
Pnigalio 
P. sp2 SCgzkd 6 f Pnig6 Kangding{] Сапа] Sichuan 2 400 Female DQ390430 
P. sp3 TBed 3 f Pnig3 Changdu[] Xizang 3 400 Female DQ390431 


For the nuclear gene ITSI[] K80 + I + G model was 
selected as the best-model by hLRTs[] - InL = 
7 457.8960[]. This model of substitution allowed for a 
bias in the rate of transitions and a substitution rate that 
varies across sites with G distribution of variable-rate 
site. Following the Bayesian analysis[] all species of 
Diglyphus clustered into one strongly supported clade 
with a high BPP[] BPP » 95[] Fig. 6[]. As with the 
analysis of CO І data[] the clade included both D. 
bimaculatus and C. ambiguus. Other species of 
Cirrospilus fell into clades A and B[] Fig. 6[]. The 
monophyly of clades A and B were strongly supported 
[] BPP » 99[]. 

For 28S D2 sequence data[] the hLRTs selected 
ће GTR + I + G as the best-fit model[] - InL = 
2 874.0681[]. Both D. bimaculatus and С. ambiguus 
clustered together with other species of Diglyphud] Fig. 
70. Cirrospilus formed 2 different clades. Diglyphus + 
C. ambiguus was monophyletic[] but the monophyly 
was not strongly supported] BPP = 620. 


4 DISCUSSION 


4.1 Phylogenetic postion of C. ambiguus and 


monophyly of the genera Cirrospilus 

Phylogenetic analyses of three gene sequences all 
supported the cluster of both D. bimaculatus and C. 
ambiguus with the other species of  Diglyphus. 
Therefore] C. ambiguus should be transferred from 
Cirrospilus to Diglyphus and referred to as Diglyphus 
ambiguus [] new combination[]. This association is 
strongly supported by the analyses of ITS1 with a high 
valu¢] BPP = 96[]. Hansson and Газа] 1996[] placed 
C. ambiguus in Cirrospilus because the species had 
the non-vaulted vertex and the metallic body coloration. 
However{] both characters were also commonly found 
among species of Diglyphus. An anterior pair of 
scutellar setae were located very close to the posterior 
margin of mesoscutum in C. ambiguus[] Zhu et al .[] 
20020] which was more typical of Diglyphus] Fig. 40 
than other species of Cirrospilus [] Fig. 3[]. 
Furthermord] the body color of species in Cirrospilus 
ranges from metallic to non-metallic. Body coloration 
did not indicate generic membership. 

Our results supported the 
Diglyphusl] but not Cirrospilus. In the Cirrospilini[] 
Cirrospilus was the most speciose genus[] Noyes[] 
2002[]. The monophyly of  Zagrammosoma was 


monophyly of 
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Majority rule 
Asp GXjxsts 2 f 1095 


C.diallus TBcd 102 f 3400 
A 
C.diallus TBcd 1 f 3400 
C.ambiguus VIETNAMhn 105 f Yr 
D.begini TBls 187 m 3650 


D.bimaculatus TBjzrl 128f 4700 %7 






D.bimaculatus TBls 130 f 3650 Y 
D.crassinervis QHdl 119 f 2900 


-C 


D.minoeus TBjzrl 101 1 M 4700 





D.pulchripes QHgem 146 f 2880 


D.isaea ОНа! 108 m 2790 


D.isaea Она! 126 f 2790 


D.pachyneurus QHdl 121 m 2900 





C.vittatus QHgem 104 f 2880 > B 
á C.variagatus SXhyhs 4 m 1475 С 
[ Eulophus sp2 TBcd 102 f 3400 
Entedon sp2 BJmtg 1280 
| Entedon sp3 SCbx 2200 


Euderus sp1 BJzhwy 101 f 100 











Fig. 5 А 5096 majority rule consensus tree based 
on CO І sequence data with Bayesian 
analysis for parasitic wasps 
Posterior probability values are shown above the branches. Letters 
indicate nodes discussed in text. Stars indicate the position of C. 
ambiguus and D. bimaculatus . The same in Figs. 6— 7. 


supported by two distinctive characters] non-metallic 
body color and a vaulted vertex. Two species[] С. 
variegatus and С.  coachellae [| Gates 2000 [T] 
supported the grouping of Cirrospilus апа 
Zagrammosoma|] as both genera share a moderately 
vaulted vertex. However[] a vaulted vertex also raises 
the possibility of a paraphyletic Cirrospilus . Ubaidillah 
et al П 2003[] analyzed 56 morphological characters of 
53 cirrospilines and showed that а monophyletic 
Cirrospilus could not be morphologically diagnosed. 
Identically[] the molecular analysis of Gauthier et al. 
[] 2000[] failed to diagnose Cirrospilus . 
4.2 Phylogenetic postion of genera Aulogymnus 
The higher taxonomic placement of Aulogymnus 
has been uncertain. When Zhu et al [] 1999[] recorded 
the genus Aulogymnus from China[] they were uncertain 
whether or not Chinese species of Aulogymnus should 
be included in the Cirrospilini because their specimens 
did not have a transverse sulcus. The genus 
Aulogymnus was placed in Ophelimini by Boucek 
[] 1988[]. Moreover[] this genus also differed from the 
other Cirrospilini in mostly having 2-4 funicular 
segmentd] Zhu et al .[] 1999[] Zhu and Huang[] 20000. 
Gauthier et а! [] 20000 revised the Eulophidae based 
on their molecular phylogeny but they did not include 
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Fig. 6 A 5096 majority rule consensus tree based 
on 1151 sequences and Bayesian 
analysis for parasitic wasps 


any species of Aulogymnus . We found that Aulogymnus 
should be included in the Cirrospilini. 
4.3 Evolution of notauli 

The Eulophidae is one of the most speciose 
families of parasitic wasps. With their diverse biology[] 
eulophids offer very good experimental model systems to 
study a variety of evolutionary questions[] Godfray[] 
1994[]. A sound phylogenetic tree facilitates an 
understanding of evolutionary history. However[] the 
high levels of homoplasy and small number of 
informative morphological characters make it difficult to 
form a hypothesis of relationships. —Fortunately[] 
molecular data cast light on the phylogeny of wasps. 
Gauthier et al.[] 2000[] proposed a phylogenetic 
framework for Eulophidae based on nucleotide sequence 
data from D2 region of 285 rDNA. 

Phylogenetic positions of D. bimaculatus and С. 
ambiguus offer an excellent opportunity to study the 
evolution of notauli patterns. The shape of notauli is not 
only an important character in eulophid taxonomy[] it is 
also associated with the ability to fly. The notauli are 
external indicators of the internal phragmata that 
separate the dorsolongitudinal and dorsoventral flight 
muscles] Michenei[] 1944[]. The phragmata form an 
internal ridge for the attachment of flight muscles 

O Wong[] 1963[]. Complete notauli indicate relatively 
more flight muscles and[] thus[] a stronger capability of 
flight. In contrast[] incomplete notauli indicate fewer 
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Fig. 7 A 50% majority rule consensus tree based on 28S sequences analyzed by Bayesian analysis 


flight muscles and a relatively weaker ability to fly. 

In the Tetrastichinae[] the notauli of all species are 
complete and extend to the hind margin of the 
mesoscutum. In contrast] the  notauli of most 
entedonines are incomplete. The notauli of euderines 
are complete and extend to the hind margin of the 
mesoscutum or sometimes disappear anteriorly. The 
phylogeny inferred from the 28S sequence shows that the 
notauli evolved from being complete to incomplete in the 
Eulophidae. The plesiomorphic condition is that of 
having complete notauli. This condition also implies 
that the capability of flight became gradually weaker as 
species of free-living eulophid ancestors specialized to 


their respective hosts . 

The complete notauli that curve to meet the apex of 
the axillae occur in five clades[] Fig. 8[]. The same 
type of notauli is also found in a few other eulophines[] 
including species of Oxycantha and Naumanniola 
ramose. Although the phylogenetic relationships among 
these species are unknown due to lack of specimens[] 
our results imply that this type of notauli evolved at 
least five times. Similarly[] incomplete notauli are found 
in four claded] Fig. 80 and[] thus] must have evolved 
independently at least four times. The multiple[] 
independent origins of derived notauli shapes suggests 
that shape of this character is phylogenetically 
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misinformative. Notauli patterns may continue to serve 
as diagnostic characteristics only when taken in concert 
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with other morphological characters . 
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Fig. 8 Evolutionary cladogram of notauli shape in the Eulophinae 


Circled notes indicate the possible formation of incomplete notauli . 


The Eulophinae is an excellent taxon for the study 
of the evolution of parasitic wasps. The phylogeny 
inferred from both molecular and morphological data can 
provide a historical context for study of evolution of 
various characters and biological features. Furthermore[] 
work in progress on the phylogeny of the Ешорћтае 
and species-level phylogenies on related genera with 
similar interspecific variation in notauli[] will greatly 
improve our understanding of the evolution of notauli[] 
particularly when more detailed biological observations 
are available for a number of species. 


Acknowledgements We thank our colleagues in the IZCAS for 


assistance with sequencing and comments on the manuscript. Dr. Yu- 














Feng Hsu[] National Taiwan Normal University[] gave us valuable 
suggestions. This work was supported by the National Natural Science 
Foundation of Chind] NSFC grant No. 30330090[] and partially by 
the National Science Fund for Fostering Talents in Basic Research 
































NSFC-J0030092[]. ZHU Chao-Dong was supported by a К.С. Wong 























Fellowship from Royal Societ] ОКО and by a scholarship from the 


Chinese Academy of Sciences to visit the British Natural History 
Museuni] Feb. 2004 – Маг. 20050. 





References 


Arias МСП Sheppard WSO 1996. Molecular phylogenetics of honey bee 
subspecies[] Apis mollifera L.[] inferred from mitochondrial DNA 
sequence. Mol. Phylogenet. Evol .[] 50 557 – 566. 

Ashmead WH[] 1904. Classification of the chalcid flies of the superfamily 
Chalcidoidea[] with descriptions of new species in the Carnegie Museum[] 
collected in South America by Herbert H. Smith. Mem. Carnegie 


SHA Zhong-Li[] et а A revised taxonomy of Cirrospilus ambiguus based on molecular systematics 
10 with notes on notauli evolution in Eulophinae 65 





Mus .0 D 400 225 – 551. 

Avise JC[] 2000. Phylogeography[] the History and Formation of Species. 
Harvard Univ. Press[] Cambridge] MA. 

Bouček Z0 1988. Australasian Chalcidoidea [] Hymenoptera [I] A 
Biosystematic Revision of Genera of Fourteen Families[] with a 
Reclassification of Species. CAB International[] Wallingford] UK. 

Buckley ТЕП Arensburger P[] Simon СП Chambers СКО 2002. Combined 
data[] Bayesian phylogenetics[] and the origin of the New Zealand cicada 
genera. Syst. Ent. 05104 – 18. 

Campbell BC[] Steffen-Campbell JD[] Werren JHO 1993. Phylogeny of the 
Nasonia species complex] Hymenoptera[] Pteromalidae[] inferred from an 
internal transcribed space] ITS2[] and 28S rDNA sequences. Insect 
Mol. Biol .[]2[] 225 - 237. 

Campbell ВО Heraty J[] Rasplus JY[] Chan КО Campbell JSO Babcock СП 
2000. Molecular systematics of the Chalcidoidea using 28S-D2 rDNA. 
0 Austin AD[] Dowton M eds. Hymenoptera[] Evolution[] Biodiversity 
and Biological Control. Australia] CSIRO. 59 – 73. 

Chen Y[] Xiao НО Fu JZ Huang DWO 2004. A molecular phylogeny of 
eurytomid wasps inferred from DNA sequences of 28S[] 18S[] 16S and 
CO І genes. Mol. Phylogenet. Evol .[ ЗШ 100 300 – 307. 

Gates MW[] 2000. A new species of Cirrospilus Westwood] Hymenopteral] 
Eulophidad] from the southwestern United States and Mexico. Proc. 
Ent. Soc. Wash .[] 102] 100 58 – 61. 

Gauthier МО LaSalle ЈО Quicke РО Godfray HCJ[] 2000. Phylogeny of 
Eulophidad] Hymenoptera[] Chalcidoidea[T] with a reclassification of 
Eulophinae and the recognition that Elasmidae are derived eulophids. 
Syst. Ent .[] 250 521 – 539. 

Gibson GAPO 1986. Evidence for monophyly and relationships of 
Chalcidoidea]] Mymaridae[] and Mymarommatidae [] Hymenopteral] 
Terebrantes[]. Can. Entomol .[] 1180 205 – 240. 

Godfray HCJ[] 1994. Parasitoids[] Behavioural and Evolutionary Ecology. 
Princeton Univ. Press[] Ottawa. 

Goldman N[] 1993. Statistical tests of models of DNA substitution. J. Mol. 
Evol .[] 360 182 - 198. 

Gordh СО 1978. Taxonomic notes on Zagrammosoma|] a key to the Nearctic 





species and descriptions of new species from Californid] Hymenopteral] 
Eulophidae[]. Proc. Ent. Soc. Wash.[] 80 300 344 – 359. 

Gordh СП 1982. Taxonomic recommendations concerning new species 
important to biological control. Int. J. Entomol .[ 10 15 – 19. 

Hall TAL] 1999. BioFdit[] a user-friendly biological sequence alignment editor 
and analysis program for Windows 95/98/NT. Nucl. Acids. Symp. 
Ser .[] 410 95 – 98. 

Hansson СП LaSalle JL] 1996. Two new eulophid parasitoidd] Hymenopteral] 
Chalcidoidea[] Eulophidae[] of Liriomyza trifoli[] Burgess([] Dipteral] 
Agromyzidae[]. Orient. Insect .[] 300 193 – 202. 

Hebert РОМ Cywinska АП Ball SI[] de Waard JRO 2003. Biological 
identifications through DNA barcodes. Proc. R. Soc. Lond. Ser. ВО 
2700 313 - 322. 

Hebert PDN[] Penton EH[] Burns JM[] Janzen DH[] Hallwachs W[] 2004. Ten 
species in one[] DNA barcoding reveals cryptic species in the Neotropical 
skipper butterfly Astraptes fulgerator. РМАЅ 100 4100 14 812 - 
14 817. 

Hoy МАО Jeyaprakash АП Morakote КО Paul KCL] Nguyen RO 2000. 
Genomic analyses of two population of  Ageniaspis  citricola 

[] Hymenoptera[] Encyrtidae[] suggest that a cryptic species may exist. 
Biol. Control] 170 1 – 10. 

Huelsenbeck JP[] Bollback JP[] 2001. Empirical and hierarchical Bayesian 
estimation of ancestral states. Syst. Biol .[] 500 351 – 366. 

Huelsenbeck JP[] Ronquist FR[] 2001. MRBAYES[] Bayesian inference of 
phylogenetic trees. Bioinformatics] 170 754 – 755. 


Ji YJO Zhang DX[] He LJ[] 2003. Evolutionary conservation and versatility of 
a new set of primers for amplifying the ribosomal internal transcribed 
spacer regions in insects and other invertebrates. Mol. Ecol. Notes[]3 

0 401 581 – 585. 

LaSalle ЈО 1989. Notes on the genus Zagrammosoma[] Hymenopteral] 
Eulophidae[] with description of a new species. Proc. Ent. Soc. 
Washington{] 910 230 – 236. 

Michener CD[] 1944. Comparative external morphology] phylogeny[] and а 

classification of the bees[] Hymenoptera[]. Bull. Am. Mus. Nat. 

Hist .[] 820 157 – 326. 

Minkenberg OPJM[] van Lenteren JC. 1986. The Leafminer Liromyza 

bryoniae and L. trifoli[] Diptera[] Agromyzidae[T] their Parasites and 

Host Plants[] A Review. Agricultural University Wageningen Papers. 

Noyes JS[] 2002. Interactive Catalogue of World Chalcidoidea[] Taxapad 

20020. Electronic Publication] CD-ROMI[J. 

Nylander ЈАП Ronquist ЕО Huelsenbeck JP[] Nieves-Aldrey Ј 2004. 
Bayesian phylogenetic analysis of combined data. Syst. Biol.[]53[] 47 
- 67. 

Peck OQ 1963. А catalogue of the Nearctic Chalcidoidea[] Insecta[] 
Hymenoptera[]. Can. Еп] Suppl.[T] 300 1 – 1 092. 

Peck O[] Boucek Z[] Hoffer A[] 1964. Keys to the Chalcidoidea of 
Czechoslovakia[] Insecta] Hymenoptera[]. Мет. Entomol. Soc. 
Can .[] 34[] 1 - 130. 

Posada D[] Crandall KA[] 1998. MODELTEST{] testing the model of DNA 
substitution. Bioinformatics[] 14] 9[T] 817 — 818. 

Ronquist ЕО 2004. Bayesian inference of character evolution. Trends Ecol. 
Evol .[] 190 475 – 481. 

Rossler YL] Debach P[] 1973. Systematic] morphology and biological 
control. Entomophaga[] 180 215 – 222. 

Sambrook J[] Fritsch EF[] Maniatis T[] 1989. Molecular Cloning[] A laboratory 
Manual. Cold Spring Harbor Lab. Press[] Plaiview[] NY. 

Schauff MEQ 1984. The Holarctic genera of Mymaridae[] Hymenopteral] 
Chalcidoidea[]. Mem. Entomol. Soc. Wash .[ 120 1 - 67. 

Schauff ME[] LaSalle J[] Coote LD[] 1997. Chapter 10. Eulophidae. 

“ Annotated Keys to the Genera of Nearctic Chalcidoidea 
П Hymenoptera”. 10 Gibson САРО Huber JT eds. NRC Research 
Press[] Ottawa. 

Simon СП Frati F[] Beckenbach АП Crespi ВО Liu НО Flook РО 1994. 

Evolution] weighting] and phylogenetic utility of mitochondrial gene 








sequences and a compilation of conserved polymerase chain reaction 
primers. Ann. Entomol. Soc. Am. 0 870 651 – 701. 

Thompson ЈОД Higgins DGO Gibson TJ[] 1994. CLUSTAL WL] improving the 
sensitivity of progressive multiple sequence alignment through sequence 
weightind] positions-specific gap penalties and weight matrix choice. 
Nuc. Acids. Res.[] 22[] 4 673 - 4 680. 

Ubaidillah RO Lasalle J[] Quicke DLJ[] 2003. Cladistic analysis of 
morphological characters in the — eulophine tribe —Cirrospilini 

П Hymenoptera[] Eulophidae[]. Entomol. Sci.[] 60 259 – 279. 

Wong HR[] 1963. The external morphology of the adults and ultimate larval 
instar of the larch sawfly[] Pristiphora erichsoni{] Htg.[[] Hymenoptera[] 
Tenthredinidae[]. Can. Е.П 950 897 – 921. 

Zhu CD[] LaSalle ЈО Huang DW[] 1999. A study on Chinese species of 
Aulogymnus Fórste£] Hymenoptera[] Eulophidae[]. Entomol. Sin.[] 60 
299 – 308. 

Zhu CD[] Huang DW[] 2002. A study of Chinese Cirrospilus Westwood 

П Hymenoptera[] Eulophidae[]. Zool. Stud .[] 10 23 - 46. 

Zhu СРП LaSalle ЈО Huang DW[] 2000. A review of the Chinese Diglyphus 

Walkef] Hymenoptera[] Eulophidae[]. Orient. Insect .[] 34[] 263 – 288 . 





udi 


66 0000 Acta Entomologica Sinica 500 





UHUUUYUUUUUUUUUUUUO 
А a ВА А 


000 70000 'О вођен W. MURPHYOOOO ™ 


01.0000000000000 10008002. 000000000 DD D 1000390 
3. Centre for Biodiversity and Conservation Biology[] Royal Ontario Museum[] 100 Queen’ s Park[] Toronto 


Ontario] M5S 260 Canada] 4. OO 00000000000000 2710180 









































































































































































































































































































































































































































































































































































































































ugguagapapapupupugggggggggggggggpapupupudmduaudu 
ummapnapaupggagpgagugnggugugnpgugugpngagagpugubnducegdg utbdubl 
Cirrospilus ambiguus L] J D. O U UU О Diglyphus bimaculatus L] D] O UO OU D] Cirrospilus D] O O O O D. О Diglyphus 
uU guBgbiudub uuu HOUUUUUUUUUUUUUUU HUUUUUYUUUUUUUUU со 
1000000 11510 080000 28502 000000 300000W000000000000000 
ugggaggagnagapapapupupugnugpnapnapnapapanpd 285020 0В600000 600 
ugggmgagagapnapnapppBppnpuaddgeourndlg udcteutl HUUUUYUUUUUUU 
HA OEN BUO A RA EE ae HEHE H EAEE 5000000 0000 0О 500 E 
ugggggaggapnppapupmugudngaggauaguu 400 

ОООоОооооооооо 0 B dgHglli DOTATA TEH HEDH 

HUUUUU Q969 HOUUUUUA 0000 0 0454-6294 2007001-0057-10 

































































